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Abstract 

The loop equation for the complex one-matrix model with a multi-cut struc- 
ture is derived and solved in the planar limit. An iterative scheme for higher 
genus contributions to the free energy and the multi-loop correlators is pre- 
sented for the two-cut model, where explicit results are given up to and in- 
cluding genus two. The double-scaling limit is analyzed and the relation to the 
one-cut solution of the hermitian and complex one-matrix model is discussed. 
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1 Introduction 



Within the topic of matrix models multi-cut solutions are of considerable 
interest. They are intimately related to the existence of multi-critical points. 
In order to reach, say, an mth critical point, ni — 1 coupling constants have 
to be introduced in the matrix potential and adjusted in the right way. This 
immediately leads to the possibility of multi-cut solutions because there can be 
as many cuts as minima of the potential. 

Unfortunately there is not much known about higher-order contributions 
in the topological expansion for multi-cut solutions. The saddle-point ap- 
proximation provides only the planar solution, whereas the full non-perturbative 
treatment with orthogonal polynomials has been successful only for one or at 
most two cuts except for special cases like degenerate minima of the potential 
Q. The reason is that the appropriate ansatz for solving the string equation 
is not known in general. The assumption for the recursion coefficients yielding 
several continuous functions in the large- AT limit does not match with the semi- 
classical analysis for higher-order potentials p. Numerical studies 0) H 
have shown the existence of instabilities in the solution of the string equation 
for a variety of different potentials. This phenomenon has been subsumed under 
the catchword of 'chaos in matrix models'. There have been attempts to explain 
the origin of these oscillations Q but a full understanding is still lacking. 

Within the framework of the third method of solving matrix models, the 
technique of loop equations Q , there has been significant progress during the 
last years. Ambj0rn et al. ^ IC, 11| have proposed a very effective scheme to 



calculate higher-genus contributions in the perturbative expansion. Making use 
of a redefinition from coupling constants to moments it allows one to determine 
all multi-loop correlators order by order in the genus expansion. 

The aim of this paper is to demonstrate how this method can be applied to 



multi-cut solutions, where the complex matrix model [12| has been considered 
for simplicity. The loop equation |^ and the starting point, the planar solution 
of the one-loop correlator, can be obtained for any number of cuts. However, 
for more than two cuts technical difficulties enter the game via a new type of 
equation determining the edges of the cuts. So the complete iterative solution 
of the two-cut complex one-matrix model presented here may be seen as a 
first step towards a possible investigation of the 'chaotic phenomena' with the 
method of loop equations. 

The paper is organized as follows. In section 2 the basic definitions are 
given. Section 3 sets up the loop equation for multi-cut correlators and its 
planar solution. In section 4 the iterative solution for two cuts is explained in 
detail, and explicit results for genus one and two are obtained. The last section 
before concluding is devoted to the double-scaling limit. 



^The approach of [ho| adopted here considerably differs from [llS 
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2 Basic Definitions 



The complex one-matrix model is defined by the partition function 

Z[N,{gi}] = e^'^ = /#td</. exp(-iV Try((/.V)) (2.1) 



with 



oo 



V(<AV) = E-('^V)' ' (2-2) 

where the integration is over complex N x N matrices. The generating func- 
tional or one-loop average is given by 

Introducing the loop insertion operator 

the generating functional can be obtained from the free energy 

W{p) = ^{p)F+-. (2.5) 
av p 

More generally one gets the n-loop correlator by iterative application of the 
loop insertion operator to F, 

W{pu . . . ,p„) = A(p„)^(p„_,) . . . ^{pi)F , n > 2 , (2.6) 

where 

, _ ^ (Tr(</>t0)fci...Tr(</>t0)fcn)^^^^^ 

W{p,,...,pn) = (2.7) 

and conn refers to the connected part. As the multi-loop correlators and the 
free energy have the same expansion 

oo -. 

W{p^,...,Pn) = 5]^M^,(pi,...,p„) , (2.8) 



9=0 



E ' (2-9) 

5=0 



relation (^) holds for each genus separately. 
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3 The Loop Equation 



In this section the loop equation for the complex matrix model [|T4| and its 
planar solution will be given. The explicit formulas are restricted to the two- 
cut case for simplicity and for consistency with the following sections. The 
multi-cut case can be found in the appendix. The origin of technical difficulties 
for more than two cuts will be also explained in this section. 

The form of the loop equation depends explicitly on the number of cuts 
of the one-loop correlator only via the contour C of the complex integral (see 
appendix^. One has 

i ^?^"'<'^> = <«'(^»' + Tf^S^ . P ^ - (3.1) 

Here, V{io) = J^JLi ^'^^"' i the support of the eigenvalue density in the two- 
cut case is o" = [—x,—y] U [y,x]. The contour C depicted in Fig.l encloses all 
eigenvalues in such a way that p can also take values on the open real interval 
between the cuts of W{p). The generalization to s cuts is obvious. 




Figure 1: the integration contour C = Ci U C2 



Now the iterative solution of ( p. ID works as follows 10, 11]. First the 
planar solution Wo{p) is determined by taking the limit iV — > 00, omitting the 
last term on the r.h.s.. This solution will then be used as a starting point for 
the iteration, which calculates Wg{p) step by step from terms of lower genera. 
For the case of two cuts, Wo{p) is given by 

with 

'^^°^('^) = , L 2 fy • (3.3) 

The extension to s cuts can be found in appendix ^ Equation (|3.2D differs 
from the one-cut solution not only by the second square root but also by a 
factor of p instead of w in the numerator. This stems from the fact that 
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depending on the number of cuts s being even or odd the complex function 
^ {lo^ — Xi) ■ ■ ■ {lo"^ — Xg) is to be defined as an even or odd function of to re- 
spectively (see also appendix P). 

From the unit normalization of the eigenvalue density it follows that 

lim W(p) = - . (3.4) 

p^oo P 

The leading asymptotic term must be accounted for already by the planar 
solution Wo{p) as 1/p does not depend on N. Imposing this condition on eq. 
(O) one finds 



1 / ^,^fcy'(^)^(o)(^) ^ A; = and 2 , (3.5) 

2 Jc 47ri 



which implicitly determines x and y as functions of the coupling constants gi. 

At this point it should be mentioned that for more than two cuts the condi- 
tion (|3.4D does not supply any more enough equations to determine all endpoints 
of the cuts. In the complex matrix model with an s-cut solution there are s such 
parameters Xi to be determined. In eq. (l3.5| ) k then runs over s, s — 2, s — 4, . . . 
down to or 1 depending on whether s is even or odd. So this yields only 
(s -|- 2)/2 or (s -|- l)/2 equations for s even or odd. The missing equations can 
be derived from the requirement that the chemical potentials between the cuts 
are equal ||l^, namely 

' dX p{X) = 0, i = l,...,s-l . (3.6) 

The eigenvalue density p(A) is given in the next section and the cfj are the 
bounded connected components of the real complement of the support a of p(A). 
Because of symmetry only the intervals cTj on the positive real line need to be 
considered, which provides the remaining (s — 2)/2 or (s — 1)/2 equations. They 
lead to a more complicated dependence of the Xi containing elliptic integrals, 
except in the case of s = 1 or 2 where they are trivially fulfilled. To see this 
the loop insertion operator is applied to eq. (|3.6|), which yields 



X /■ V^^^ """i) ■ ••• • [a seven 

= Ml'>—^{p) / dX^ -2 2^ <^ 3.7 

^ dV Ja-j (A2 -xf) [ ^ ^ ^^"^ • 

The moments m[^^ are to be defined in the next section. Together with ^(p) 

of eq. ( |3.5[ ) the relations ( |3.6D determine the derivatives ^(p) for i = 1, . . . , s, 
which are needed in the iterative process. However, their complicated structure 
makes it hard to see whether a scheme for calculating higher genera can still be 
established. 

Coming back to the iterative solution of the loop equation it turns out that 
after the insertion of the genus expansion eq. (|2.8D , Wg (p) is determined by the 
following equation 



9-1 



d 



{k-2Wo{p))Wgip) = J2Wg,{p)Wg-g,(j)) + ^ip)Wg.i(j)) , g>l . (3.8) 
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Here, the linear operator K is defined by 

mp) - f ^ (3.9) 

JC 4:111 P'^ — LO^ 

Wg{p) is now expressed completely in terms of Wg,{p) with gi < g on the r.h.s. 
of eq. ( |3.8| ). Hence the next step is the inversion of the operator {K — 2Wo{p)) 
acting on it. In contrast to the one-cut case this operation will involve zero 
modes contributing to Wg{p) which have to be fixed. 



4 The iterative solution 
4.1 Change of variables 



In analogy to ||Il|] it is convenient to change variables from the coupling 
constants gi to the moments and in the following way 

^ r du^r{u^)fH^ ^ ^^^^ 

The advantage of these new variables is that, for given genus, Fg and the multi- 
loop correlators Wg{pi, . . . ,Pn) depend only on a finite number of moments 
instead of the infinite set of couplings. Moreover, the mth multi-critical point 
can be characterized by the vanishing of the first m — 1 moments or J^. 
This happens whenever extra zeros of the eigenvalue density 

p{X) = -\ImWo{X)\ = -3-|M(A)|7(x2-A2)(A2-y2) , xea. (4.2) 
vr 47r * 

occur at either x"^ or y^. The analytic part M(A), given in appendix must 
develop these extra zeros. Using eqs. ( [B.9D and ( |4.1| ) it can be seen that the 
moments defined above provide an expansion of M(A), i.e. 

oo 

M(A) = ^(MfcA(A2-x2)'^'i + JfcA(A2-yy-i) . (4.3) 
fc=i 

Calculating the moments in terms of the coupling constants yields 



Ml 



k 



9k+i + {{k + ^)x^ + ^y^^ gk+2 + • • • , 



Jk = 9k+i +({k + ^)y^ + ) gk+2 + ■■■ ■ (4.4) 
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4.2 Inversion of {K - 2Wo{p)) 

In order to proceed in solving equation ( |3.8D it is necessary to find a set of 
basis functions for the operator {K — 2Wo{p)) acting on Wg{p) as follows 

{K - 2Wo{p)) x^-Hp) = J^^^n>l, 

iK-2Woip))i^^''\p) = j^^^,n>l. (4.5) 

As the main result, Wg{p) will then be expressed in terms of these functions, 

3g-l 

Wg{p) = 5:4")x^"^(p) (4.6) 

n=l 

where the coefficients Ag^^ and Bg"'^ only depend on the moments and 
and on and y'^. 

Note that eq. ( |4.5| ) does not yet completely determine x^^\p) 'ip^"'\p) 
due to a non-trivial kernel of {K — 2Wo{p)). Because of eq. (|3.4D only terms 
asymptotically of order p~^, k > 2, can contribute to Wg{p) for g > 1. In the 
one-cut case this made the definition of the basis unique in a simple way. The 
argument excluded that the zero mode Xj \Jp^ — of {K — 2Wo(p))one-c«tj 
which is asymptotically of order i, could be added to Wg{j)). However, here 
this is no longer the case as 

Ker(i^-2Wo(p)) = Span |p(^(°)(p), (4-7) 

and 

~ — . (4.8) 

p p" 

Hence such a term can be added to Wg(p) in every step of the iteration. How 
shall y}^"'\p^ and ip^^^p) be fixed such that Wg{p) is uniquely determined? By 
definition (see eq. (|2.5| )) Wg{p) can be written as a total derivative 

Wgip) = -^{p)Fg, g>l. (4.9) 

In order to satisfy this equation, the p-dependence of x^^\p) ^-iid ip^'^\p) must 
be completely absorbed in terms of the type ^(p), ^(p)) ^f^(p) ^(p)- 
Consequently the basis functions x^'^Kp) 'il)'^"'\p) must be linear combina- 
tions of these, the coefficients of course again depending on moments, and 
y^. In this manner the zero mode contributions to the basis, which are indeed 
necessary, become uniquely fixed. The final result, proven in appendix]^, takes 
the following form 



i^^'^Hp) ^ ( (4"^(p) - E i^^'\p)Jn-k) - E i^^'\p)Jn-k+i ) ,(4.10) 



k=l 
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where 



^^\P) - ^^n^n>l. (4.11) 



4.3 The iterative procedure determining Wg{p) 

In order to solve eq. ( |3.8| ), one first has to compute ^(p)VFo(p) = Wo{p,p)- 
To do so it is convenient to rewrite the loop insertion operator 

^<^' = i^^f^' + ^f^'il^ + ^wlj. (4.12) 

where 

("3) 



For evaluating ^(p) and ^(p), eq. ( [4. 121) is applied to eq. (|3^) which yields 
Now with eq. ( p^ ^(p)VFo(p) 

can be calculated. Using the relation 



the definitions ( |4.ip and deforming the contour of the remaining integral to 
infinity, one gets 

W0{P,P) = 16(^2 _ ^2)2(^2 _ ,2)2 ' (4-16) 

where 

d = x^-/. (4.17) 

Before one can make use of the basis functions eq. ( ^OD for {K - 2Wo{p)) one 
has to decompose the r.h.s. of eq. j ) into fractions of the form (p^ — x^)^"^ 
and {p^ — y^)~" . Doing so the coefficients Ag"'' and in ([4.6| ) can now be 
identified for genus 1. 



^ 16 8 d ' 2 16 

5f ) = - + -4 , = - (4.18) 

1 16 8 d ' 2 16 ^ ' 



This simple relationship with the basis functions is spoiled for more than two cuts. 
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It is clear how to carry on the iteration process. However, in order to calculate 



'^'{p)Wg{p) it is convenient to rewrite the loop insertion operator once again, 



-iP) = y'-^{p)^ + y'-^iv)^ + 'Ap)— + 'Ap)— 

dV^^' ^ dV ^^'dMn ^,dV ^^'dJn dV ^^'dx^ ^ dV ^^'dy^ ' 

n=l n=l " ^ 



(4.19) 



-^ip) = -in + -^)<^^-^'Hp)--^T.yi^U'Hp)-M,^ 



+(n + i)M„+i^(p) . (4.20) 

With the definition ([4.1| ) of the moments and the form of ^ (p) in eq. ( [4.12| ) the 
relation is easily verified. The corresponding result for ^^{p) can be obtained 
by interchanging x"^ <-> and ^ J^. Looking back at the original loop 
equation (|3.1|) it is obvious that the free energy and the multi-loop correlators 
for all genera should be invariant under such an interchangement . 

Using the loop insertion operator in the form of eq. ( |4.19| ) a lengthy calcu- 
lation yields the following result for g = 2: 

(1) ^ 5 y^Ms 5 xV3 1 x^MaJa 1 M2J2 

2 32d3(Mi)3 32 d3(Ji)3 64(i3(Mi)2(Ji)2 + 128(i2(Mi)2(Ji)2 

49 y'^{M2f ^ 49 x^(J2)^ ^ 11 x'^M2 11 y'^J2 



256d3(Mi)4 256(i3(Ji)4 128d4(Mi)2Ji l2Sd^Mi{Ji 
67 y2j\^2 67 j;2j2 3 M2 3 J2 



128d4(Mi)3 128d4(Ji)3 128d3(Mi)2Ji 128d3Mi(Ji)2 
201 201 x2 57 57 1 



256d5(Mi)2 256 d5(Ji)2 64d5MiJi 128^4 Mi Ji ' 

(2) ^ 5 j/^Ms 1 x^M2J2 49 y''{M2f 
2 32d2(Mi)3 + 128(i2(Mi)2(Ji)2 256 d2(Mi)4 

3 x2M2 ^ 1 2;2J2 3 J2 67 y2^2 



128d3(Mi)2Ji 16d3Mi(Ji)2 128d2Mi(Ji)2 128d3(Mi)2 
57 2;2 201 y2 ^5 ^ 
128^4 Ml Ji ~ 256d4(Mi)2 128^3 Mi Ji ' 

(3) _ 5 x^Ms 49 x'^{M2f 5 XV2 67 x2M2 



2 32d(Mi)3 ^ 256 d(Mi)4 128 d2Mi( Ji)2 ^ 128 d2(Mi)3 

49 M2 15 x2 201 x2 189 1 



(4) _ ^ - I 



128d(Mi)3 128(i3MiJi 256d3(Mi)2 256 (i2(Mi)2 
49 x2M2 189 x2 105 1 



2 128d(Mi)3 256 d2(Mi)2 ' 256d(Mi)2 

(5) _ 105 x2 



2 256d(Mi)2 



5^'^ = A'^''{Mk^ Jk,x' ^y') , i = l,...,5 . (4.21) 
W2(p) is then obtained by inserting eqs. ( [4.10| ) and ( 4.21 ) into eq. (^). 
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4.4 The iterative procedure for Fg 



As it was already pointed out in section 4.2 and proven in appendix |C| the 



basis functions x^"^^{p) c^iid 'ip^"'\p) can be written as linear combinations of 
derivatives with respect to V{p). This requirement was imposed in order to fix 
the basis in a unique way. For genus 1 the explicit result reads 



X 



y2 dV 

which can easily be verified from the definitions. The relation ( [4. 91) then allows 



to calculate Fg for any given Wg{p). So from eq. (4.22) in combination with 



eqs. (4.18) and ( [4. 61) Fi can now be read off up to a constant 



Fi = In(Mi) - 1 ln( Ji) - ^ ln(<i) . (4.23) 

Continuing in the same manner and rewriting the basis as being sketched in 
appendix the genus two result is obtained using eq. ( [4.21| ) after some tedious 
work, 

35 M4 35 J4 43 M3 43 J3 



384(i(Mi)3 384d(Ji)3 imd'^{Mif 192(i2(j^)3 
29 M3M2 29 J3J2 1 M2J2 



128(i(Mi)4 128d(Ji)4 64d2(Mi)2(Ji 

11 t 1^ 



21 (M2)3 21 (J2)3 11 (M2)2 11 (J2^2 



+ 



160d(Mi)5 160d(Ji)5 40d2(Mi)4 40d2(Ji)4 
181 M2 181 J2 3 M2 3 J2 

~l~ TTTTT 7Q / T \ Q TTT 7q / 71 \ o t ~I~ 



480(i3(Mi)3 480(i3(Ji)3 64d3(Mi)2Ji 64d3Mi(Ji)2 
181 1 181 1 5 1 

~480d4(Mi)2 ~ 480d4(Ji)2 ~ led^MiJi ' ^' ^ 

From eq. ( [4.6D and the procedure described just above it should have become 
obvious that Fg depends on at most 2{3g — 2) moments. 

Another remarkable fact is that the calculated Fi and F2 exactly coincide 
with those of the one-cut hermitian matrix model described in |11] when the 



identification of the moments and Jk and the difference d is made (using 
the same notation). This coincidence away from the double-scaling limit cannot 
be merely pure coincidence. However, since the loop insertion operators are 
distinct the loop correlators will differ as it can already be seen for Wi{p) and 
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W2{p)- Relating the two models it has been mentioned in ||T^ that the complex 
matrix model corresponds to a hermitian matrix model with a general potential 
where the eigenvalues are restricted to be positive ^ It seems that the repulsion 
at the origin (A = 0) is reflected only in the multi-loop correlators but not seen 
in the free energy. 

Finally a comparison to the one-cut solution of the complex matrix model 



presented in |10] can be made by letting y — > 0. The moments defined here and 



there then translate into each other. 

Mk{gk+i ^ 9k) Ik 

Jk ^ Mk (4.25) 

As one might have expected, the results for the free energy and the multi-loop 
correlators do not match except for the universal 2-loop correlator Wo{p,p). For 
a given set of couplings the phase boundary between the one- and two-cut phase 
can be formulated. The free energy or the correlators can be compared on the 
boundary, inspecting the order of the phase transition and critical exponents. 



4.5 The general structure of Fg and Wg{p) 

The main result for Fg from the iterative solution of the loop equation can 
be written in the following way 

Fg = («i..-«A.;/?i...AI«,A7)c; — rf7(Mi)"(Ji)/^ g>2. 

(4.26) 

Here the brackets denote rational numbers and a, (5 and 7 are non-negative 
integers. The summation-indices Oj and [3j take values in the interval [2, 3(7 — 2]. 
For every genus g, Fg contains only finitely many terms with a finite number 
of moments. In particular Fg depends on at most 2(3g — 2) different moments. 



This structure can either be proven along the same lines like in [10| or becomes 



clear when performing the first steps of the iteration. In perfect analogy to the 



one-cut case of the hermitian matrix model [11| several relations between the 
indices and powers in eq. ( [4.26| ) can be derived. 

First of all because of the symmetry y^, Mk Jk the following holds 

(ai . . . afc;/3i . . . Aja, /3, 7)^ = {-If {(3i .../?;; ai •• • ak\f3,a,-f)g . (4.27) 

Defining 

Nm = k-a , Nj = l-P , (4.28) 

it is true that 

Nm < , Nj < . (4.29) 



^This can be seen after diagonalisation when the eigenvalues of (p^cf) are considered to be 
A > instead of with A real here. 
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The invariance of the partition function Z =exp(^g A^^ ^^-^3) under the rescal- 
ing N kN and p{\) \p{^) for each genus yields 

Nm + Nj = 2-2g . (4.30) 

The rescaling k'^N and gj — > k^^^gj implies 

Mj k^-^Mj , Jj kp-^Jj , 

because of eq. ( [4.4[ ). For Fg this reads 

k I 

5^(a,-l)+^(/?, -l)+7 = 45-4 . (4.32) 
i=i j=i 

In the double-scaling limit in the next chapter further relations of this type will 
be derived allowing to decide which terms in eq. ( ^.26| ) will survive. 

Turning to Wg{p) the coefficients Ag^'^ and Bg""^ have a similar structure 



like in eq. ( 4.261 ) since Wg{p) follows from Fg by applying the loop insertion 



operator. In particular, 

(4.33) 

for g > 1 and analogously for Bg . The only difference is that the indices a, 
and f3j lie in the interval [2,3g — n]. Because of the same genus expansion eq. 
(4.28) is also valid here. So with the definition ( [4.1[1| ) of the basis functions 



Wg{p) depends on at most 2(3(7 ~ 1) nioments. 

5 The double-scaling limit 

In the conventional double-scaling limit all matrix models belonging to the 
same universality class should be equivalent. Consequently all differences origi- 
nating from the multi-cut structure should vanish in this limit. Having explicit 
results at hand for the one-cut hermitian and the one- and two-cut complex 
matrix model ([^,jl^]), this can be checked as an example. 

In the following the double-scaling limit will be performed for only; the 
one for is then easily obtained. For the mth multi-critical model the couplings 
are adjusted such that 

x"^ = xj.— aA™ , 

= x1 + avr (5-1) 

and 7/^ does not scale. Here a is the scaling parameter, which becomes zero at 



the critical point. As it had been already mentioned in section |4.1| the eigenvalue 
density then develops m — 1 extra zeros at and hence 



Mk ~ a""-^ , k = l,...,m-l , (5.2) 
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whereas the J^'s do not scale. The resulting contribution to the free energy has 
the well known behavior 

Fg - a(2-29)(-+|) ^ ^ > 1. (5.3) 
Making use of eq. (|]2|) for Fg leads to 

k 

^(m - Ui) - a{m - 1) > m{2 -2g) - g + l . (5.4) 
1=1 

Since in the scaling limit the free energy should look the same for all multi- 
critical models it follows that 

Nm > 2-2g, 

k 

^ 3g-3. (5.5) 

1=1 



From eqs. (l430|) and (|I29| ) the equality sign holds in the first equation. Only 



terms for which this is true also in the second line will contribute in the double- 
scaling limit. For example the genus two contribution to F will then be 

pid.s.i.) ^ 35 M4 29 M3M2 21 {M2f 
2 384d(Mi)3 128d(Mi)4 160 d{Mi f ' ^ ' ' 

Switching to the one-loop correlator Wg{p) the behavior of the basis func- 
tions in eq. ( ^.1C1| ) also has to be analyzed. Equations ( |5.1| ) and ( |5.2| ) lead 
to 

V^(")(^,A) ~ . (5.7) 

The known genus g contribution 

Wg{Tr,A) - a(i-2s)(m+i)-i (5.3) 

(n) 



together with eq. ( 4.33| ) requires for the Ag terms that 



k 11 

^(m — ttj) — a{m — 1) — m — n -\ — > m(l — 2g) — g . (5-9) 

i=l 

The same argument as above then yields 

Nm > 2-2g, 

k 

- 1) ^ 3<7 - n - 1 , (5.10) 

i=l 

where again only terms obeying equality appear in the scaling limit. A similar 
computation for Bg^-' reveals Nm > 1 — 2g, which clearly cannot be fulfilled as 
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an equation together with eq. ( [4.30| ) because Nj < 0. So all the Bg"'^ terms will 
disappear in the double-scaling limit. Finally the non- vanishing coefficients of 
W2{p) are given as an example. 

^(3) _ 49 x2(M2)2 5 xlMs 



A. 



2 

(4) 

2 128 4(Mi)3 



256 4(Mi)4 32dc{Miy 
49 x^Ms 



4'^ = — (5.11) 

In W2(p) the x^-dependence cancels out because in the basis the second sum 
is suppressed in the scaling limit. This reproduces exactly the result for the 
one-cut hermitian matrix model in |l^], where the equivalence to the one-cut 
complex matrix model had already been proven. 



6 Conclusions 

It has been shown how the powerful method of iteratively solving the loop 
equation by Ambj0rn et al. ||9|, |l^, 11| generalizes to the complex matrix model 



with more than one cut present. The loop equation for an arbitrary number 
of cuts was derived and solved for the one- loop correlator in the planar limit. 
In principle the procedure to find the genus g contribution is clear also for the 
multi-cut type. Nevertheless, for more than two cuts a new kind of equation 
determining the edges of the cuts enters and renders the computation technically 
much more involved. 

The iterative scheme was then explicitly presented for the two-cut model, 
and results for genus one and two were obtained away from the double-scaling 
limit. Relations to the one-cut solution of the hermitian and complex matrix 
model were discussed, in particular in the case of the double-scaling limit and 
when the two cuts merge. 

In order to attack the problem of instabilities in multi-cut solutions termed 
'chaos in matrix models' a more complicated cut structure has to be examined, 
e.g. the hermitian model with two cuts for an arbitrary potential or simply 
just with three or more cuts. Up to now the instabilities have only been found 
within the picture of orthogonal polynomials. The hope is that they can be 
found also in the framework of loop equations and that a deeper understanding 
especially concerning correlation functions can be obtained in this way. These 
open problems are left subject to further investigations. 



Acknowledgements: I would like to thank P. Adamietz, J. Ambj0rn, O. Lech- 
tenfeld and J. Plefka for helpful discussions. 
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A Derivation of the loop equation for the s-cut so- 
lution 



The matrix integral eq. ( p.lD is invariant under the foUowing transformation 



The functional determinant and the change of the action is then given by 

^ (^1 + 2ep(TV-2^^)2^ , (A.2) 



y(<AV) - V{ct>U) + le J'^^^^ V'icpU) , (A.3) 

where 

oo 

V'i^U) EE 5:5n(</>V)"-' • (A.4) 
n=l 

The invariance of Z then reads 

(Tr^^)') - ^ (Tr^^F'(0t0)^ = , (A.5) 



— (p^ (p / \ 



which is already almost the loop equation using eqs. (2^) and (2.7) 



Now the density of the eigenvalues can be formally introduced as 

PNW^^(T.S{^->^i)) ■ (A.7) 



\i=l 



With a = Uf^^Uj being the support of our s-cut solution the explicit dependence 
on the number of cuts s enters the l.h.s. of eq. (|A.6|). 



Tc 4:TTi uj^ — — uj^ 

duj 2uPV^^) 

— uP- 



f :r-W{uj) \' . (A. 
Je 4:711 — 



Here C encloses all cuts without containing ztp, which generalizes the two-cut 
case depicted in Figure |l] of section ^. With the following change of notation 
eq. (|3.1D is finally obtained. 



oo 



2ujV'{uj'^) = 2 ^ gni^^""^ = V'{uj) . (A.9) 

n=l 
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B The planar solution of the loop equation 

In the limit of A*" — > oo the loop equation (|3.lD becomes 
Jc 4:111 p'^ — U)^ 

Deforming the contour C to infinity and using the fact that W{j)) and V (j)) are 
odd functions by definition one gets the following contributions from the poles 
at ibp and oo 

{W,{p)f = \V\V)W,{V) ^ i ^^^P4^oM. (B.2) 
2 J Coo 4vrz — 

The solution of this quadratic equation for Wo(p) of course reads 



Q^v) = i 7-:T±im^) (B.4) 



Wo(j>) = \v'{p) ± yi(y'(p))2 , (B.3) 

with 

duj ujV'{uj) _ 
fCoa 47ri — op- 
to be calculated for any given potential with finitely many couplings. Making 
an ansatz for a solution with s cuts Wq{j)) looks the following 

W^o(p) = \ {y\v)-M{p)^{p^-x'\)...{p'-x'^S) , (B.5) 
where M{p) is an analytic function. So 

M{p) = <^(°)(p)(F'(p)- Wo(p)) , (B.6) 

remembering 

y/(p2 - xf) . . . (p2 _ xj) V ^'V V 

(B.7) 

Now for s even (odd) c^^^^ (p) is defined as an even (odd) complex function of p 
and consequently 



Mip) = UM{p) + {-iy+'M{-p)) 



2 

1 / duj J 1 (-l)^+^\ ^ 
= - j —M{^)[ +^-^—] ■ (B.8) 

Reinserting eq. ( [B.6| ) into the integral in eq. (|B.§| ) the term proportional to 
Wo(p) drops out because of its asymptotic eq. (|3.4D and hence 

X / duj 2V'(uj) jf]\, . \ p seven ,^ 
/Coo 47r« u;2 - p2 ^ 1 a; ^ odd . 



Plugging this into eq. ( p.5[ ) again after a similar calculation Wq{p) can then be 
expressed as 



1 / (hj V'{uj) 



= 2 fc iS^^^M^ • I ^ ^ odT. (^-1°) 
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C Determination of the basis 



As being described in section 4.2 our aim is to find a basis which can be 
expressed completely in terms of moments, and as well as total derivatives 
^(p) of them. The results for the latter were already given in eqs. (3) and 
( 4.20| ). It will be shown by induction that the basis defined in eq. ( 4.1C| ) can be 
expressed in the conjectured way. The starting point was made in eq. (|4.22| ). 
Now assume that this holds for all x^^\p) ■• A; = l,...,n — 1. Rearranging eq. 
( [4.10 ) shows that the same is true then for the (l)x\p) 



k-l 



X 



1. (C.l) 



1=1 



1=1 



It follows with eq. ( 4.20| ) that the remaining term in x^"'\p) can also be rewrit- 
ten in the desired way. 



dMn- 



n 



dV 



-ip) 



^n-l (_^)n+fc-l 
2 ^ 



{4'\p)-mAp) 



k=l 



d 



•n-k 



(C.2) 



The proof for is going exactly along the same lines. The fixing of the 

basis is unique because the zero mode alone cannot be written as a derivative 
with respect to V{p). 

What remains to show is that x^^\p) ^"^^ il)^"'\p) really form a basis like in 
eq. (4.5). First define in analogy with the one-cut case 



x^'^Kp) 



1 



1 



n-l 



(C.3) 



which does not contain the zero modes yet. It is easily proven by induction 
that 

{k - 2Wo{p))x^-\p) = _^^2)n ' " > 1 ' 
holds. Now the zero modes are added to x^"'\p)i 



n—1 



Ml \p^ a;-" p 



„2n 



J2x^''Hp)Mn-k+i] , (C.5) 



so eq. ( |C.4| ) without tilde is still valid. Finally the equivalence of eq. ( p. 51 ) to 
the form in eq. ( 4.10 ) is again shown by induction. Proceeding in the same 
way for ip^'^^p) completes the proof of this section. 
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